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Swinging of red blood cells under shear flow 
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We reveal that under moderate shear stress (777 m 0.1 Pa) red blood cells present an oscillation of 
their inclination (swinging) superimposed to the long-observed steady tanktreading (TT) motion. A 
model based on a fluid ellipsoid surrounded by a visco-elastic membrane initially unstrained (shape 
memory) predicts all observed features of the motion: an increase of both swinging amplitude and 
period (1/2 the TT period) upon decreasing 777, a 777-triggered transition towards a narrow 777- 
range intermittent regime of successive swinging and tumbling, and a pure tumbling motion at 
lower 777- values. 
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A human red blood cell (RBC) is a biconcave flattened 
disk, essentially made of a Newtonian hemoglobin solu- 
tion encapsulated by a fluid and incompressible lipid bi- 
layer, underlined by a thin elastic cytoskeleton (spectrin 
network) [l[ . The complex structure of RBCs and their 
response to a viscous shear flow have a great influence on 
flow and mass transport in the microcirculation in both 
health and disease Q ■ The full understanding of this re- 
sponse requires a direct comprehensive observation of cell 
motion and deformation, and a model for deducing the 
cell intrinsic properties from its behavior in shear flow, ft 
is generally admitted that the two possible RBC move- 
ments are the unsteady tumbling solid-like motion 
and the drop-like 'tanktreading' motion for higher shear 
stresses, where the cell maintains a steady orientation, 
while the membrane rotates about the internal fluid, as 
reported respectively for RBCs suspended in plasma or in 
high- viscosity media and submitted to high shear stresses 
@) EL HI El- However, the RBC movement at smaller 
shear rate and close to the tumbling-tanktreading tran- 
sition, has not been fully explored. Moreover, the actual 
state of deformation of the elastic skeleton either in the 
flowing or in the resting RBC is still conjectural ("shape 
memory" problem) [3] ■ Most models [a, Q derive from 
the analytical framework of Keller and Skalak (KS) [3], 
which treats the RBC as a fluid ellipsoidal membrane 
enclosing a viscous liquid. Although this model qualita- 
tively retrieves the two modes of motion, it does not cap- 
ture the observed shear-rate dependency of the tumbling- 
tanktreading transition. In particular, the model does 
not account for the possible elastic energy storage in- 
duced by the local deformations of the cytoskeleton dur- 
ing tanktreading. Approaches including membrane elas- 
ticity are either restricted to spherical resting shapes be- 
cause of analytical complexities [lo| . or propose encour- 
aging but still limited numerical analysis on tanktreading 
elastic biconcave capsules 



Here, we reveal a new regime of motion for RBCs under 
small shear flow, characterized by an elastic capsule-like 
oscillation of the cell inclination superimposed to tank- 
treading that we name swinging. We develop a model, 
which predicts both swinging and the shear-stress depen- 
dency of the tumbling-tanktreading transition. It demon- 
strates the existence of the elastic shape memory in the 
membrane. 

Direct measurements of cell orientation with respect 
to the flow direction (angle 9) and cell shape (lengths of 
the long and small axis of the cell cross-section, ai and 
a2 respectively) are provided from side-view microscopic 
imaging in a vertical plane parallel to the plane of shear 
[13]. We varied the wall shear rate 7 (in the range 0- 
5 s" 1 ) and the outer viscosity r\ by suspending RBCs 
in various solutions of dextran (concentration 6%, 7.5% 
or 9% w/w and viscosity 22, 31 and 47 mPa.s respec- 
tively). Therefore, the wall shear stress is varied in a 
range from to 0.25 Pa. For the highest values of the 
external shear stress 770-7, tanktreading is observed. It is 
characterized by i) a quasi-stationary cell shape with in- 
significant deformation (maximum variation of ai < 5%), 
ii) rotation of the membrane, revealed from the motion 
of small carboxylated beads stuck to the membrane (Fig. 
[23) and iii) an oscillation of the cell inclination about 
a mean value ranging from 6° to 25° (Figs. UJP and [2j 
see [lH) at a frequency equal to twice the tanktreading 
frequency (Fig. QJ3). Such characteristics are not seen 
on tanktrea ding vis cous lipid vesicles [13, 14 1, neither 



predicted [3, O, 0, E3] ■ This oscillation is however ob- 
served for non-perfectly spherical elastic millimeter-scale 
capsules [3, EH (Fig. 03) or protein-coated drops [2fj| . 
and in numerical simulations on biconcave elastic shells 



111 ]. We explain this phenomenon by assuming RBC 
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shape memory. We state that the local elements of the 
composite membrane (cytoskeleton and lipid bilayer), in- 
cluding the elements which form the rim and the dim- 
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FIG. 1: Units [r?i] = [?7 ] = [?7m]=mPa.s, [/i]=Pa and .[7]=s . (A) Schematic drawing of a tanktreading ellipsoid in a shear flow. 
(B) Rotation of a bead (diameter 1/im) stuck on the membrane of a RBC with (7 = 6, T)o = 47). Time sequence of Is. (C) 
RBC swinging : (7 = 1.33, rj — 47). Time sequence of 2s. 
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FIG. 2: Same units as in Figure 1. (A) Orientation versus 
the normalized time 7t for various cells from top to bottom: 
(7 = 1.8, Vo = 22), (7 = 2.6, ^ = 31), (7 = 6.6, Vo = 47). (B) 
Orientation versus -yt for: (o) a RBC with (7 = 0.8, r\ — 47); 
solid line from (Eq. EJ with (r/ m = 1120, fi m = 0.42); (□) a 
polymeric capsule from JTiJ with (7 = 18, r\ = 964). Solid 
line from (Eqns. [5]) with the surface moduli [r\ m .e = 0.085 
mPa.s/m, /i m .e = 0.675 mPa.m), and (ai = 278.8/xm, 02 = 
az — 170.8 /im) obtained from the size at rest Ro = 224.8 /im 
of the capsule and its mean deformation during flow: D — 
(01 - a 2 )/(a 1 + aa ) « 0.12 at 7 = 18. 



pies are not equivalent and are not strained in the bi- 
concave resting shape. They do consequently not store 
elastic energy. Thus, during tanktreading, the elements 
which form the rim at rest rotate about the stationary cell 
shape to reach the dimples after rotation and reciprocally. 
These elements are then locally strained and store elastic 
energy. Both local deformation and energy storage are 
periodic: each time the elements of the membrane make 
a 7r-rotation, they retrieve their initial shape and are no 
more strained. We emphasize that the periodic storage of 
energy requires a non spherical unstrained state for the 
RBC. Otherwise, the membrane elements will tanktread 
without modifying the global state of stress of the cell, 
preserving the steady nature of the tanktreading move- 
ment. 



In order to derive tractable equations of motion, we use 
the KS model. We consider an oblate ellipsoid filled with 
a viscous liquid and delimited by a viscoelastic 3D thin 
membrane, which includes the lipid bilayer and the un- 
derlying cytoskeleton 2l|. The membrane elements are 
prescribed to rotate along elliptical trajectories parallel 
to the shear plane, with a linear velocity field 22[ given 
by: V\ = -uj(ai/a 2 )x2,v 2 = u>{a 2 /ai)xi,v 3 = 0, where 
w and u> are the phase angle of a membrane element and 
its instantaneous frequency of tanktreading respectively 
(Fig. The KS equation for RBC motion is obtained 
by stating that at equilibrium, the total moment exerted 
by the external fluid on the cell vanishes (First equation 
in EqsO below). In addition, the movement satisfies the 
conservation of energy, i.e. the rate of dissipation of en- 
ergy in the cell must equal the rate at which work is 
done by the external fluid on the cell. KS calculated 
both rates assuming viscous energy dissipation in the 
cell. We add to this latter contribution the elastic power 
stored in the periodic elastic strain of the cytoskeleton 
24j : P e i — J n Ti(a : D)dO, where f2 is the membrane 
volume, D the eulerian strain rate tensor derived from 
the KS velocity field and a the shear stress tensor in the 
membrane; a is computed from the local deformation of 
the membrane due to tanktreading, assuming a simple 
Kelvin- Voigt viscoelastic material: a = 2r/ m T) + 2/x m E, 
where E is the Euler-Almansi strain tensor obtained from 
the KS velocity field. After some algebra, P e i writes as 
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where r\ m and \i m are the membrane viscosity and the 
shear modulus respectively. Conservation of energy pro- 
vides a constraint on the allowable RBC motion and 
yiel ds a second differential equation (for more details see 
[131). The two coupled equations are: 
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where 9 is the time derivative of the cell inclination, 
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fx, h an< i h are geometrical constants and V is the 
RBC volume (same definition as in [9]). The limiting 
case /i m = r] m = corresponds to KS. We numerically 
solved the equations using the following set of parame- 
ters for RBCs : ai = a3 = 4/im, &2 = 1.5/im, f2 = E.e 
[iBj], where £ is the oblate ellipsoid area and e=50 nm 
is the membrane thickness j26[. r\i is fixed at the physi- 
ological value of 10 mPa.s [6| and rj m is adjusted in the 
range 0.7-2 Pa.s [6j. We obtain 0(i), 0(i), w{t) and u(t). 
Suitable couples of /x m and r\ m values were found to re- 
produce exp erimental RBC oscillations as seen in Fig. 
[2(3 (see [13(). One example of experimental capsule os- 
cillation extracted from [19J is presented as well in Fig. 
[2)3. An insight of experimental and numerical swinging 
curves is provided from three parameters: the magnitude 
A9 = 9 max — 8 m i n , the mean angle 9 and the period T osc 
(or the frequency f osc ) of oscillation (Fig. [2)3) . The 7- 
variations of these parameters are illustrated for one red 
blood cell in Fig. [3K together with a numerical solution 
of Eqns [2] While 9 decreases for decreasing 7 down to 
0, T osc and A9 increase until 7 reaches a critical value 
7~ below which the cell tumbles at least once. Besides 
direct observations (Fig. QJ3), we illustrate the factor 2 
of proportionality between f osc and Co which relates the 
movement of swinging of a cell to the movement of tank- 
treading of its membrane, by reporting on a same graph 
in Fig. [5J3, variations of / osc /2 and u> versus 7, for differ- 
ent swinging RBCs observed in our experiment and for 
different tanktreading RBCs observed in the litterature 
at higher 7 [!, Q . Fig. [3j3 shows indeed the continuity of 
the two regimes even though measured on different type 
of movement. Therefore, all experimental characteristics 
are well captured by the model. Indeed, by treating the 
elastic contribution as a small perturbation in the second 
equation of (2), valid in the linear part, one recovers the 
steady KS solution at the order of the parameter in 
front of s!n(2w), while at the first order one finds that 9 
oscillates at twice the tanktreading frequency (linear in 
7) and A9 scales as ([i m /j). 

Contrarily to KS prediction, we observe that the tran- 
sition of movement from tumbling to tanktreading (re- 
spectively tanktreading to tumbling) is induced by tuning 
up (respectively down) the applied shear rate. This tran- 
sition is illustrated in Fig. (4j\. Its more striking feature, 
predicted and experimentally observed is the existence 
of a regime of movement where the cells present succes- 
sively swings and tumbles at a given 7 (Fig. |4}3). The 
model gives the 7-range ([7^,7^]) where this regime ex- 
ists. Given the experimental constraints, it is not easy to 
follow the cells sufficiently long to observe a large series of 
tumbles and swings. We define over a time scale of ~ 20 
s, the shear rate corresponding to a change in movement 
from swinging to tumbling (7,5) with decreasing 7, and 
from tumbling to swinging (7?) with increasing 7. The 
difference 7?" -7,5, that we call hysteresis (for instance 
for two different RBCs at rj =31 mPa.s : 7< = 0.47 
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FIG. 3: Same units as in Fig. 1. (A) Experimental data on 
a single RBC at r) = 22: (•) AO, (■) T osc , (□) 1 tumbling 
period value and (o) 6 versus 7. Corresponding curves of the 

model with y m = 0.38 and rj m = 700: ( — ) AO, ( ) T osc 

and ( ) 9. (B) f osc /2 versus 7: (.) rj = 22; (♦) r/o = 47; 

(■) r) = 31; (o) Single RBC, r) = 22. (x) Tanktreading 
frequencies tl) taken from p] with rj = 35, (+) from [5| with 
rj = 18, 31, and 59, (E2) from [| with r/ = 70. The line is a 
guide for the eyes. 
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FIG. 4: Same units as in Fig. 1. (A) The transition from 
swinging to tumbling induced by decreasing 7 is associated 
with a transient localized deformation (r] =A7, 7=2.66). Time 
sequence of Is. (B) (o) Successive swinging and tumbling at 

(j)o = 22, 7=1.526); ( ) numerical calculus with (r) = 22, 

7=1.526, (im = 0.454, r\ m = 700). (C) Theoretical shear- 
stresses of transition versus fi m with rjm = 1000: (o) r] "/c\ 
(x) ?7 7+ and (a) Vo(jc~ ~ > + ). 

s- 1 , t> = 1 s" 1 and j< = 0.8 s _1 , j> = 1.73 s" 1 
respectively) has the same order of magnitude than the 
theoretical 7-range [jj ,7^] of intermittency. Finally as 
it is seen by requiring the second term in the second equa- 
tion of (2) to be of the same order of magnitude as the 
first part, the critical shear rate should scale as fJ-m/Vo- 
It is indeed numerically observed (Fig. HP). Both criti- 
cal values of r\oit an d loic are mainly governed by the 
RBC elastic contribution for given cell geometry and may 
provide an average determination of fj, m by observing a 
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rheology measurements. 
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FIG. 5: Distribution function of critical shear stresses of tran- 
sition for increasing ( A / > ) and decreasing (7 < ) shear rates. 



large sample of RBCs (Fig. HP). Distribution functions 
of r/oj^ and i r\oif' measured on a large RBC sampling 
are shown in Fig. [5] They characterize the natural vari- 
ability of the RBC elastic modulus and illustrates the 
additional hysteresis effect. From Fig. [5]and Fig. HP, we 
find that fi m ranges in the interval 0.14-2 Pa. By setting 
the 2D shear modulus fi m .s = Mm- e we obtain values 
ranging from 0.07 to lxlO -7 N/m (e = 50 nm) below 
that usually reported 0, [l(| ■ We also find a comparable 
difference on /J, m ,s with that reported for the elastic cap- 
sules by Walter et al This underestimation likely 
originates from the major simplifications we made in or- 
der to obtain simple analytical equations allowing the full 
understanding of the physics of the problem: i) simplis- 
tic constitutive equations, ii) approximate KS velocity 
field, which may overestimate membrane deformations. 
In particular, Tran son Tay et al Q suggested that the 
Secomb-Skalak area conserving velocity field [23| would 
lead to a 70%-increase of the membrane viscosity com- 
pared to that derived from the KS model, iii) treat- 
ment of deformations from a 3D description of the mem- 
brane of RBCs and capsules although these systems form 
mainly 2D-shells 2l|. However, the main interest of this 
tractable model is to understand the role of the various 
physical parameters on the movement. For example, for 
given external viscosity and shear rate, A9 is not much 
sensitive to values of r\i and r\ m taken in the physiolog- 
ical range. A9 is essentially fixed by /i m , and thus the 
measurement of the amplitude of swinging as a function 
of the shear stress may provide a complementary method 
to accurately determine fj, m on single flowing RBCs in a 
given sample. 

In conclusion, the swinging movement and the shear- 
stress triggered transition of motion of RBCs demon- 
strate the existence of their shape memory and is a sig- 
nature of their membrane shear elasticity. Despite its 
simplicity, our model provides a good description of the 
observed behaviors and we believe that a more sophis- 
ticated model should allow an easy and sensitive deter- 
mination of individual RBC mechanical properties. Fi- 
nally, such experimental approach coupled to the refined 
model might be applied to a wide variety of soft shells 
, 1 1 QL l20j ] and holds promises for applications in surface 
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